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At optical frequencies the radiation produced by a source, such as a laser, a black 
body or a single photon source, is frequently characterized by analyzing the temporal 
correlations of emitted photons using single photon counters. At microwave frequen- 
cies, however, there are no efficient single photon counters yet. Instead, well developed 
linear amplifiers allow for efficient measurement of the amplitude of an electromagnetic 
field. Here, we demonstrate how the properties of a microwave single photon source 
can be characterized using correlation measurements of the emitted radiation with such 
detectors. We also demonstrate the cooling of a thermal field stored in a cavity, an 
effect which we detect using a cross-correlation measurement of the radiation emitted 
at the two ends of the cavity. 



In quantum optics [T], the single photon detector is a 
versatile tool to explore the properties of radiation emit- 
ted from a variety of classical and quantum sources. At 
optical frequencies these detectors easily produce a 'click' 
when a single photon impinges on them. Recording the 
statistics of such events, one is able to measure not only 
the average of the number of emitted photons n = (aJa) 
but also higher order statistical correlations between the 
emitted photons. Here, a' and a are the creation and an- 
nihilation operators of the radiation detected in a given 
mode. At lower frequencies, such as in the microwave 
frequency domain, however, no efficient single photon 
detectors exist to date as photons at these frequencies 
carry an energy that is orders of magnitude less than 
at optical frequencies. Instead, at microwave frequencies 
linear amplifiers followed by an instrument recording a 
voltage, such as an oscilloscope, are commonly used to 
detect small amplitude electromagnetic fields. Such am- 
plifiers have several orders of magnitude gain but also 
add noise to the signal. Nevertheless, valuable informa- 
tion about the properties of quantum radiation sources in 
the microwave domain can be acquired. With sufficient 
averaging, electric field amplitudes, which are propor- 
tional to the sum and difference of the field operators a 
and cv, can be detected in this way even on the single 
photon level [2J. We demonstrate here that by recording 
the time series of the detected signal - instead of its time 
average - and using efficient digital signal processing, we 
can extract characteristic correlation functions of radia- 
tion generated by microwave frequency emitters such as 
a single photon source or a black body. 

For our experiments we have realized a microwave fre- 
quency single photon source in a superconducting elec- 
tronic circuit, similar to the one presented in [5J. In 
our circuit, we coherently and controllably couple a sin- 
gle qubit to a high quality resonator to create an indi- 
vidual photon on demand. The superconducting trans- 



mon qubit [3] used in this experiment is characterized by 
its maximum Josephson energy Ej tinax w 14.4 GHz, its 
charging energy Eq ~ 500 MHz and its energy relaxation 
and dephasing times in excess of a few hundred nanosec- 
onds. The transition frequency of the qubit is flux tun- 
able using both a quasi-static magnetic field generated 
with a miniature coil and an on chip transmission line 
to generate nanosecond time scale flux pulses. By inte- 
grating our qubit into a superconducting coplanar trans- 
mission line resonator of frequency v T = 6.433 GHz and 
quality factor Q = 2060 we couple it strongly to a sin- 
gle mode a of the radiation field stored in the resonator. 
This approach is known as circuit quantum electrody- 
namics [3] and allows to study in exquisite detail the 
interaction of quantum two-level systems with quantized 
radiation fields. 

Applying a phase controlled truncated gaussian mi- 
crowave pulse of variable amplitude A r and total dura- 
tion t-c — 18 ns to the qubit biased at a transition fre- 
quency of P q — 6.933 GHz, we prepare a superposition 
state \tp q ) = a\g) +/3|e) between the qubit ground \g) and 
excited states |e). The superposition is characterized by 
the two complex probability amplitudes a = cos(# r /2) 
and j3 — sin(# r /2)e l * which are parameterized by the 
polar (Rabi) angle 6 T and the phase angle </>. We char- 
acterize the prepared qubit state using a pulsed disper- 
sive measurement of the resonator transmission [5] and 
clearly observe Rabi oscillations in the qubit population 
P e versus the amplitude A r (Fig. [iji). After the qubit 
state preparation, we apply a current pulse of controlled 
amplitude and duration to the flux bias line to tune the 
qubit transition frequency into resonance with the res- 
onator frequency z/ r . We time-resolve the resonant vac- 
uum Rabi oscillations of the coupled system at a fre- 
quency of g/n — 118 MHz by dispersively measuring the 
qubit state after it has been tuned back to the frequency 
is q strongly detuned from the resonator (FigfTV>). Adjust- 
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FIG. 1: Control of qubit state and qubit/resonator in- 
teraction, a, Qubit excited state population P c versus Rabi 
pulse amplitude A r . The Rabi rotation angle 9 r is indicated 
on the top axis. The used pulse shape is shown in the inset, b, 
P a versus qubit-resonator interaction time t VT with the used 
flux pulse shape in inset. Note that the interaction time can 
experimentally not be made shorter than 1.9 ns due to finite 
rise time of the pulses. Dots are data, lines fits to theory. 



ing the qubit-resonator interaction time to half a vac- 
uum Rabi period i vr = n/g — 4.2 ns, we coherently map 
the qubit state |^ q ) to an equivalent superposition state 
|V> C ) = ot\Q) + f3\l) of the |0) and |1) photon Fock states 
stored in the resonator mode a. Similar techniques have 
been used to prepare and measure a wide range of intra- 
cavity photon superposition states in recent experiments 
both with superconducting circuits [6] and with Rydberg 
atoms [7]. 

In a next step, we characterize zero and one photon 
superposition states by measuring the quadrature ampli- 
tudes of the microwave fields emitted from both ends of 
the cavity, see schematic in Fig. [2j For this purpose we 
have realized a symmetric circuit QED setup in which the 
output fields are detected independently and simultane- 
ously at both ports of the cavity. It is interesting to note 
that this setup is essentially equivalent to one in which 
the output of a radiation source is investigated at the two 
output ports of a beam splitter [8] ■ Our scheme comprises 
two independent detection chains similar to the one pi- 
oneered by Gabelli et al. [5]. Each chain consists of a 
cold amplifier with gain Gb, c ~ 33 dB and noise temper- 
ature lN(b,c) ~ 4.5 K followed by a two stage heterodyne 
detector in which the signal is down converted from the 
resonator frequency to 25 MHz in an analog stage and to 
d.c. in a digital homodyne stage. This allows for the mea- 
surement of the electric field at both ends of the resonator 
characterizing the radiation emitted into the modes la- 
beled b and c. In this way, we extract the complex enve- 
lope St c {t) of the amplified electric field E^J described 
by 



(1) 



where the real and imaginary part of St, jC (i) are the two 
field quadratures in the frame rotating at the resonator 
frequency and Nb, c (t) is the noise added by the ampli- 
fier. Using input-output theory |10j . one can show that 
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FIG. 2: Schematic of the experimental setup. Resonator 
in which a two level system emits a photon into the mode a. 
Linear amplifiers with gain Gt, c couple symmetrically to the 
cavity outputs amplifying the radiation in output modes b 
and c. Homodyne detectors (HD) extract both quadratures 
of the output fields and feed results to field programmable 
gate array (FPGA) based digital correlation electronics. 



the full information about the intra-cavity mode a can be 
extracted from a measurement of the propagating modes 
b and c [8]. As a first example, we present a measure- 
ment of the time dependence of both quadrature ampli- 
tudes of the electric field of the output mode b at one 
end of the resonator. This gives us access to the ex- 
pectation value of the annihilation operator of the cavity 
field (Sb(t)) oc (a(t)) [8]. Similar measurements were pre- 
sented in Ref. [2J, where the cavity photon was created 
by Purcell limited spontaneous emission. Figure|3^, shows 
the real part of (S b (t)) versus time t after the prepara- 
tion of the photon superposition state \ip c ) characterized 
by the qubit Rabi angle 6 r used for its preparation. We 
find excellent agreement with the expected average field 
quadrature amplitude (a) oc sin (9 r )/2 (Fig. |3f;), where 
we find the largest signals for the superposition states 
l^c + ) = (|0> + |1»/V2 and |C) = (|0> - |1»/V§ pre- 
pared using 6 r — it/ 2 and 37r/2, respectively. As ex- 
pected from the uncertainty principle, the Fock states 
|0) and |1) prepared with 9 r = and tt, respectively, 
do not show any quadrature amplitude signals (Fig. [3^) 
since their phase is completely uncertain. For all of the 
above measurements, the overall global phase of the sig- 
nals is adjusted such that the imaginary part of Sb(t) 
is equal to zero which therefore is not displayed. We 
also note that the amplifier noise averages to zero in the 
quadrature amplitude measurement. Moreover, the time 
dependence of all measurement traces is well understood. 
For the state |?/>+) (Fig. |3p), for example, the character- 
istic decay time is given by twice the cavity decay time 
2r K = Q/nv r = 102 ns. The rise time, which should ide- 
ally be t vr = 4.2 ns, is limited by the bandwidth 14 MHz 
of our detection scheme. 

In our measurement scheme, we simultaneously record 
the time dependent quadrature amplitudes Sb, c (t) de- 
tected at both ports of the cavity for each single photon 
that we generate. This is realized using a two channel 
analog-to-digital converter (ADC) with a time resolution 
of 10 ns. Based on the measurement record of each event, 
we can then calculate any expectation value, such as av- 
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FIG. 3: Quadrature amplitude and cross-power mea- 
surements, a, Measured time dependence of cavity field 
quadrature amplitude of mode b for zero and one photon su- 
perposition states characterized by the Rabi angle 8 T (left 
axis) and generated state (right axis), b, Single quadrature 
trace at 9 T = in a corresponding to (|0) + 1))/ a/2- c, 
Dependence of maximum quadrature amplitude on 9 r . d, 
Measured time dependence of cross-power between modes b 
and c for the same preparation as in a. e, Single power trace 
at 8 T = 7T in c corresponding to |1). f, Dependence of maxi- 
mum cross-power on 9 r . Blue dots are data, red line is theory, 
see text. 



erages, products or correlations, that can be expressed in 
terms of the detected output signals Sb, c (t)- Processing 
this data in real time using field programmable gate array 
(FPGA) electronics allows us to efficiently extract infor- 
mation even in the presence of substantial noise added 
by the amplifier. 

As a first experiment taking advantage of this scheme, 
we have measured the expectation value of the instanta- 
neous power (S£(t)Sb(t)} emitted into the output mode 
b with the cavity mode a prepared in the Fock state |1) 
(not shown). It is important to note that we have digi- 
tally calculated the power by multiplying the quadrature 
amplitudes Sb(t) for each single photon pulse instead of 
using a diode as a power meter in which the detection 
and the averaging is realized within the detector [2J. As 
we will see in the following, this digital procedure is very 
versatile as it also allows to calculate cross-powers and 
cross-correlations digitally. In the direct power measure- 
ment, the detected noise power of the amplifier dominates 
by a factor of about 2000 over the single photon power 
which is still observed using sufficient averaging (similar 



as in Fig. [3)3). From the background noise we determine 
the system noise temperature T^^s) w 17 K of our de- 
tection chain with respect to the output of the resonator. 
^N(sys) is substantially higher than the noise temperature 
of the amplifiers because of absorption in the cables and 
insertion loss of components in the detection chain. 

Calculating the cross-power (Sl(t)S c (t)} between the 
two output modes instead of the direct power emitted 
into just a single output mode, we can reject most of the 
noise added by the amplifiers, demonstrating the versa- 
tility of our digital scheme. The detected cross-power is 
related to the average photon number in the cavity as [5] 



P(Ni 
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where P(Nb c ) is the power of correlated noise between 
channels b and c. In these measurements, the detected 
noise cross-power is a factor of 1/500 smaller than the 
direct noise power of each amplifier as the two detec- 
tion chains add predominantly uncorrelated noise. The 
residual correlations are of technical origin, such as in- 
sufficient isolation of the two amplifier chains, correlated 
digitizer noise and residual resonator thermal noise due 
to incomplete thermalization of the resonator inputs. We 
have characterized the measured cross-power of our sin- 
gle photon source for the same set of cavity superposi- 
tion states as used for the quadrature amplitude mea- 
surements (Fig. [3b) . We find excellent agreement of the 
temporal evolution of the cavity photon number in de- 
pendence on the preparation angle of the photon state 
(ala) oc (1 — cos# r )/2 (Fig. [3^) . The maximum cross- 
power is measured for the Fock state |1) (8 r = it) and the 
minimum power for the |0) state (9 r = or 2tt) (Fig.[3[l). 

Finally, we have characterized our single photon source 
using time-dependent first-order cross-correlation mea- 
surements of the two output modes of the resonator 



r (1) (r)= I (St(t)S c (t + r))dt. 



(3) 



For this purpose, we generate a train of 40 single pho- 
ton pulses, each created using the procedure described 
above, with a pulse separation of t p — 512 ns which is 
much greater than the qubit and cavity decay times. To 
remove the background, we subtract the measured corre- 
lation function rii (r) in the resonator steady-state from 
the signal acquired when performing the photon state 
preparation sequence rW(r). From the recorded quadra- 
ture amplitude data, we calculate in real-time 



r«(r)-rW(r)cxG«(T). 



(4) 



which gives us access to the first-order correlation func- 
tion GW(t) = / (a^(t)a(t+T))dt of the resonator field 0. 
To measure each trace in Fig. [4^,, 128 x 10 6 trains of 40 
photons were prepared in a specific state and G^(t) was 
calculated in real time using our FPGA based electron- 
ics, corresponding to more than 1 terabyte of data that 
have been evaluated in approximately 1 hour. 
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FIG. 4: Correlation function measurements, a, Time 
dependence of first-order correlation function of the cavity 
field G"' for indicated states. Data is offset for improved 
visibility, b, Correlation function at r = and r = nt p 
versus 9 r . Dots are data, lines are theoretical predictions, 
c, Enlarged trace for (9 r = from a displaying cavity field 
cooling. 



The correlation function data G^(r) (Figj4p,) is char- 
acterized by a set of peaks that are separated by the 
repetition time t p of the single photon source. The am- 
plitude of G^ at r = and r = nt p , representing the 
correlation between a pulse i and i + n, depends in a 
characteristic fashion on 6 r . For the Fock state |1) (at 
9 r = 7r), the correlation function G^(0) is at a max- 
imum and vanishes at G^^ntp) as there is no coher- 
ence between photons emitted from the source at differ- 
ent times. In fact, G^(0) oc (a^a) oc |/3| 2 = sin((9 r /2) 2 
oscillates sinusoidally with the preparation angle, as it es- 
sentially measures the average photon number of the gen- 
erated field (Fig. [4Jd) . For photon superposition states, 
the expectation values of (<v) and (a) of subsequently 
generated photon states have non-vanishing values, as 
discussed before. Since photons from different repetitions 
of the experiments are uncorr elated, G (1) (n* P ) 

cx (a 1 ') (a) 

which has a finite value and oscillates at half the period. 
Thus, G^(nt p ) oc \a/3\ 2 — sin(# r ) 2 /4 is maximized for 
the states \ip+) and |$r) (Fig-gja). 

The observed features of the first-order correlation 
function are characteristic for the non-classical radia- 
tion generated by our deterministic single photon source. 
Considering the well controlled procedure implemented 
for generating the single photon pulses any other pulsed 
coherent or thermal field can essentially be ruled out 
as being the source of the measured correlations. This 
indicates for our scheme that the quantum mechanical 
properties of the radiation are not lost in the detection 
chain but are fully retained in the form of statistical in- 
formation of the data and can be extracted by appropri- 



ate single-shot analysis. Measurements of the second or- 
der correlation function G^ would provide unambiguous 
proof of the quantum character of the field, independent 
of any prior knowledge about its source pQ. So far, the 
long integrations times required due to the presence of 
noise added by the amplifiers have prevented their suc- 
cessful measurement. 

Interestingly, our time-dependent first-order correla- 
tion function measurements allow us to observe cooling of 
the resonator field through its interaction with the qubit 
in its ground state. If the thermal occupation of the qubit 
excited state is substantially smaller than the thermal oc- 
cupation of the resonator, the qubit can absorb a photon 
from the resonator during the interaction time t vr . The 
qubit then emits the photon into the environment at a 
different frequency thereby cooling the resonator field. 
The cooling is observable as a pronounced dip in the mea- 
sured G^(t) around r = 0, see Fig. [4J;. An analysis of 
the size of the dip in the correlation function allows us 
to extract the thermal background population of the res- 
onator ?ibg = 0.07 corresponding to a field temperature 
of Tb g = 115 mK which is in agreement with indepen- 
dent measurements of 7b g from the vacuum Rabi mode 
splitting spectrum |llj . Analyzing in detail the time de- 
pendent correlation function measurements of the single 
photon source in the presence of a weak thermal back- 
ground field, we find excellent agreement between our 
data and theory [5] for all prepared field states (see solid 
lines in Fig. 

Our experiments clearly demonstrate that correlation 
function measurements based on quadrature amplitude 
measurements are a powerful tool to characterize quan- 
tum properties of propagating microwave frequency ra- 
diation fields. Even in the presence of noise added by 
the amplifier, efficient data processing techniques allow 
for the measurements of higher statistical moments of the 
fields. When better, possibly quantum limited, amplifiers 
|12) become available the demonstrated techniques may 
help to enable the full tomography of propagating radi- 
ation fields. Furthermore, the flexibility of circuit design 
and the high level of control achievable in circuit QED 
will enable a variety of future experiments with quantum 
microwave fields for basic research and applications. 

While preparing this manuscript, we became aware of 
related work by Menzel et al. [15] . 
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